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We prove a general rigorous lower bound for VF(A, q) = exp(S'o(A, q)/kB), 
the exponent of the ground state entropy of the g-state Potts antiferromagnet, 
on an arbitrary Archimedean lattice A. We calculate large-(7 series expan- 
sions for the exact Wr{A,q) = q^^W(A,q) and compare these with our lower 
bounds on this function on the various Archimedean lattices. It is shown that 
the lower bounds coincide with a number of terms in the large-g expansions 
and hence serve not just as bounds but also as very good approximations to 
the respective exact functions ^^(A,^) for large q on the various lattices A. 
Plots of Wr{A,q) are given, and the general dependence on lattice coordina- 
tion number is noted. Lower bounds and series are also presented for the 
duals of Archimedean lattices. As part of the study, the chromatic number 
is determined for all Archimedean lattices and their duals. Finally, we report 
calculations of chromatic zeros for several lattices; these provide further sup- 
port for our earlier conjecture that a sufficient condition for Wr{A,q) to be 
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analytic at 1/g = is that A is a regular lattice. 
05.20.-y, 64.60.C, 75.10.H 



I. INTRODUCTION 



Nonzero ground state disorder and associated entropy, 5*0 7^ 0, is an important subject in 
statistical mechanics. One physical example is provided by ice, for which the residual molar 
entropy is ^0 = 0.82±0.05 cal/(K-mole), i.e., Sq/R = 0.41±0.03, where R = NAvoghs 1- 1- 
Indeed, residual entropy at low temperatures has been observed in a number of substances, 
including nitrous oxide, NO, carbon monoxide, CO, and FCIO3 (a comprehensive review 
is given in Ref. [^). In these examples, the entropy occurs without frustration, i.e., the 
configurational energy can be minimized. In magnetic systems, two examples are provided 
by the Ising antiferromagnet on the triangular and kagome lattices here, the ground 

state entropy does involve frustration. A particularly simple model exhibiting ground state 
entropy without the complication of frustration is the g-state Potts antiferromagnet (AF) [0 
on a lattice A, for q > x(A), where x(^) denotes the minimum number of colors necessary 
to color the vertices of the lattice such that no two adjacent vertices have the same color. 
As is already evident from the foregoing, this model also has a deep connection with graph 
theory in mathematics |]T^, since the zero-temperature partition function of the above- 



mentioned g-state Potts antiferromagnet on a lattice A satisfies Z{A, q,T = 0) paf = P{-^, q), 
where P{G, q) is the chromatic polynomial expressing the number of ways of coloring the 
vertices of a graph G with q colors such that no two adjacent vertices (connected by a bond 
of the graph) have the same color; hence, the ground state entropy per site is given by 
So/ks = \nW{A,q), where W{A,q), the ground state degeneracy per site, is 

l^(A,g)=JimP(A„,g)^/" (1.1) 

Here, A„ denotes an n- vertex lattice of type A (square, triangular, etc.), with appropriate 
(e.g., free) boundary conditions. Given the above connection, it is convenient to express 
our bounds on the ground state entropy in terms of its exponent, M^(A, q). Since nontrivial 
exact solutions for this function are known in only a very few cases (square lattice for g = 3 
triangular lattice |T^, and kagome lattice for g = 3 []T3|), it is important to exploit and 



extend general approximate methods that can be applied to all cases. Such methods include 
rigorous upper and lower bounds, large-g series expansions, and Monte Carlo measurements. 
Recently, we studied the ground state entropy in antiferromagnetic Potts models on various 
lattices and obtained further results with these three methods []TB|- pD|. 



In the present paper we achieve a substantial generalization of our previous studies. 
Among other things, we obtain a general rigorous lower bound on the (exponent of the) 
ground state entropy that applies for all Archimedean lattices (for definitions, see below) 
and find further examples of lattices where this lower bound coincides to many orders with a 
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large-g series expansion of the respective W function. This agreement is particularly striking 
since, a priori, a lower bound on a function need not coincide with any, let alone many, of the 
terms in the series expansion of the function about a given point. We also present calculations 
of chromatic zeros for a number of lattices and show that they support a conjecture that we 
have made earlier [^. The reader is referred to Refs. [l^- for further background and 
references. 

Although the full set of Archimedean lattices is not as much discussed in the physics 
literature as the subset of three homopolygonal (i.e., monohedral) ones, viz., square, trian- 
gular, and honeycomb, other Archimedean lattices are of both theoretical and experimental 
interest. For example, the kagome lattice is of current interest because of the experimental 
observation of compounds whose behavior can be modelled by quantum Heisenberg an- 
tiferromagnets on this lattice, including SrCr8_xGa4+a;Oi9, (SCGO(x)) and, recently, 
deuteronium jarosite, (D30)Fe3(S04)2(OD)6 [1^. Related to this, one of the reasons for 



interest in the kagome lattice is that on this lattice the quantum Heisenberg antiferromagnet 
has a disordered ground state with finite entropy and frustration [|23|. Indeed, it has been 
known for a number of years that although the Ising antiferromagnet exhibits ground state 
entropy, frustration, and zero long range order on both the triangular and kagome lattices, 
the greater degree of disorder on the kagome lattice is evidenced by the fact that while the 
correlation length diverges as T ^ on the triangular lattice p4|, it remains finite on the 



kagome lattice |25 



II. SOME GRAPH THEORY BACKGROUND 

In this section we include some basic definitions and results in graph theory that will be 
necessary for our work. A graph (with no loops or multiple bonds) is defined as a collection 
of vertices and edges (bonds) joining various vertices. In strict mathematical terminology, a 
graph involves a finite number of vertices; the regular (infinite) lattices which we will consider 
here are thus limits of graphs, with some appropriate (e.g. free) boundary conditions. The 



chromatic polynomial P(G, g), defined above, was first introduced by Birkhoff p6[ and has 



been the subject of intensive mathematical study since [27|- |28|, [^- |12|. The zeros of the 



chromatic polynomial P{G, q) are denoted the chromatic zeros of G. A graph G that can be 
colored with q colors, i.e., has P{G, q) > 0, is said to be g-colorable. The chromatic number 
x{G) of a graph G is the minimum (integer) q such that one can color the graph subject to 
the condition that no two adjacent vertices have the same color, i.e. such that P{G, q) > 0. 
The graph G is then said to be x(G')-chromatic. The graph is uniquely x(G)-chromatic 
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if and only if P{G,q = x{G)) = ql- A uniquely g-cliromatic graph is g-partite. These 
definitions can be extended to the regular infinite lattices of interest here if one defines the 
latter as the limit as the number of vertices n — oo of finite-n lattice graphs with appropriate 
boundary conditions, such as free boundary conditions. For bipartite lattices one can also use 
periodic boundary conditions that preserve the bipartite property. Similarly, for a tripartite 
lattice, one can use periodic boundary conditions that preserve the tripartite property, etc. 
Henceforth, when we refer to a finite lattice of type A, it is understood that appropriate 
boundary conditions are specified. 

An Archimedean lattice is defined as a uniform tiling of the plane by regular polygons 
in which all vertices are equivalent [^. Such a lattice is specified by the ordered sequence 



of polygons that one traverses in making a complete circuit around a vertex in a given 
(say counterclockwise) direction. This is incorporated in the mathematical notation for an 
Archimedean lattice A: 

A = iUpr) (2.1) 

i 

where in the above circuit, the notation p^" indicates that the regular polygon pi occurs 
contiguously times; it can also occur noncontiguously. Because the starting point is 
irrelevant, the symbol is invariant under cyclic permutations. For later purposes, when a 
polygon Pi occurs several times in a non-contiguous manner in the product, we shall denote 
Qi^s as the sum of the Oj's over all of the occurrences of the given in the product. There 
are eleven Archimedean lattices: 

{A} = {(3^), (4^), (6^), (3^ ■ 6), (3=^ ■ 4^), (3^ ■ 4 ■ 3 ■ 4), (3 ■ 4 ■ 6 ■ 4), (2.2) 



(3-6 ■3-6), (3-122), (4.g.i2)^ (4-82)} (2.3) 

Of these lattices, three are homopolygonal, also called monohedral, i.e., they only involve one 
type of regular polygon: (3^) (triangular), (4^) (square), and (6^) (hexagonal or honeycomb). 
The other eight are heteropolygonal, i.e., involve tilings with more than one type of regular 
polygon. The kagome lattice is (3 ■ 6 ■ 3 ■ 6). Some illustrative pictures of the heteropolygonal 
lattices are helpful for understanding the results of the present paper. Accordingly, we show 
in Figs. the (3^ ■ 4 ■ 3 ■ 4), (3 ■ 6 ■ 3 ■ 6) (kagome), (3 ■ 12^), and (4 ■ 8^) lattices. Pictures 
of others are given in Ref . p9[ ; see also Ref . ||30| . 
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FIG. 1. Section of the (3^ • 4 • 3 • 4) Archimedean lattice with paths and C!^ used in the proof 
of our lower bound indicated with darker lines. See text for discussion. 
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FIG. 2. Section of the (3 • 6 • 3 • 6) Archimedean lattice. The point A labels a vertex that is 
referred to later in the text. 
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FIG. 3. Section of the (3 • 12^) Archimedean lattice. 
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FIG. 4. Section of the (4 • 8^) Archimedean lattice. 



The geometric constraint that the internal angles of each polygon adjacent to a given 
vertex sum to 27r is 



The solutions to this equation for integral ai^s and pi, together with the requirement that a 
local patch be extendable to tile the plane, yield the Archimedean lattices. The degree A 
of a vertex of a graph G is the number of edges (bonds) that connect to this vertex. For a 
regular (infinite) lattice, this is the same as the coordination number. For an Archimedean 
lattice ( |2.1| ), the coordination number is 




(2.4) 



(2.5) 
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Of course, for a finite lattice with free boundary conditions, the vertices on the boundary have 
lower values of A than those in the interior; this will not be important for our rigorous bounds, 
which pertain to the thermodynamic limit on an infinite lattice. For a homopolygonal lattice 
A = {p"'), the constraint ( p.4|) relates the coordination number to p according to 

A = a = , for A = (p") (2.6) 

p — 2 

which can also be written in the symmetric form A^^ + p^^ = 1/2. The girth g{G) of a 
graph G is the length of a minimum circuit on the graph. Hence, for an Archimedean lattice 
(p. ID, g = mm{pj}. The number of polygons of type Pi per site is given by 

_ Np. per V _ dj^s 7\ 

^" ~ iV ~ V ^ ' 

^^v per Pi fi 

The set of homopolygonal Archimedean lattices is invariant under the duality transformation, 
which interchanges 0-cells (vertices) and 2-cells (faces) and thus maps {p"") —>■ (a^)- When 
one applies the duality transformation to the other eight Archimedean lattices, the resultant 
lattices are not Archimedean. 



The dual of the Archimedean lattice, often called a Laves lattice [p9| , pl| , is defined by 
listing the ordered sequence of vertex types specified by their degrees, Vi = Aj, along the 
boundary of any polygon: 



ALaves= [U^-'] (2.8) 



where in the above product, the notation v^* indicates that the vertex Vi with degree f j = Aj 
occurs contiguously 6j times; it can also occur noncontiguously. As with Archimedean lat- 
tices, because the starting point is irrelevant, the symbol is invariant under cyclic permu- 
tations. When a vertex of type Vi occurs several times in a non-contiguous manner in the 
product, we shall denote bi^s as the sum of the foj's over all of the occurrences of the given Vi 
in the product. There are eleven dual Archimedean lattices. Just as all vertices are equiva- 
lent on an Archimedean lattice, all faces are equivalent on its dual; these are comprised of 
a single type of polygon, p, which, however, does not in general have sides of equal length. 
Note that 

P = T.Ks (2.9) 

i 

and the girth (/([Hi^r^']) = P- The dual of the Archimedean lattice (HiPD is [Oj^r^'] with 
Vi = Pi and hi = Qi. Among Laves lattices, three involve only vertices of one type and are 
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denoted homovertitial and are equivalent to the respective three homopolygonal Archimedean 
lattices; these are [3^] = (6'^) (honeycomb), [4"^] = (4^) (square), and [6^] = (3®) (triangular); 
the other eight are heterovertitial, i.e. involve vertices of at least two types. An example of 
a heterovertitial Laves lattice is [4 ■ 8^], the union-jack lattice, shown in Fig. ^ 



FIG. 5. Section of the [4 • 8^] dual Archimedean (Laves) lattice. 



For these duals of Archimedean lattices, the formula for the number of p-gons per site is 

b- 

(2.10) 



Although different vertices on a heterovertitial Laves lattice have different degrees, it will 
be useful for later purposes to introduce the notion of an effective degree or coordination 
number, defined by Ag// = limy^oo where V and E denote the vertices and edges of the 
lattice graph. This is given by 
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Ae// = UpP 

This completes our brief review of necessary definitions from graph theory. 



(2.11) 



III. CHROMATIC NUMBERS FOR ARCHIMEDEAN AND DUAL 
ARCHIMEDEAN LATTICES 

As part of our work on ground state entropy, it has been useful to calculate the chromatic 
numbers x(^) fo^^ the Archimedean lattices and their duals, and to determine whether or 
not these lattices are uniquely g-chromatic for q = x{^)- As background, we recall that from 
the four-color theorem by duality, it follows that xiG) < 4 for any planar graph |P2|| . 
We find that for all of the nonbipartite Archimedean and dual Archimedean lattices, x = 3 
except for the [3 ■ 12^] Laves lattice, for which x = 4. The latter result follows from the fact 
that this lattice contains K4 subgraphs, and xi^m) = (Here, Km is called the complete 
graph on m vertices, defined such that each vertex is joined by bonds to every other vertex.) 
We summarize our results in Tables H and ||, together with some properties of these 



lattices that will be relevant to our calculations of bounds and series. 
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TABLE I. Chromatic number x(-^) determination of uniqueness or non-uniqueness of col- 
oring if g = x(A) for the Archimedean lattices. UQC stands for "uniquely g-chromatic" , Y,N for 
yes and no. 



A 


X(A) 


UQC 


(36) 


3 


Y 


(4^) 


2 


Y 


(63) 


2 


Y 


• 6) 


3 


Y 


(33 • 42 ) 


3 


N 


(32 • 4 • 3 • 4) 


3 


N 


(3 • 4 • 6 • 4) 


3 


N 


(3 • 6 • 3 • 6) 


3 


N 


(3 • 122) 


3 


N 


(4 • 6 • 12) 


2 


Y 


(4 • 82) 


2 


Y 
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TABLE II. Chromatic number x{^) ^i^d determination of uniqueness or non-uniqueness of col- 
oring if g = x(A) for the heterovertitial duals of Archimedean lattices. UQC stands for "uniquely 
qi-chromatic" . 



A 


X(A) 


UQC 


[63] 


3 


Y 




2 


Y 




2 


Y 


• 6] 


3 


N 




3 


N 


[32 • 4 • 3 • 4] 


3 


N 


[3 • 4 • 6 • 4] 


2 


Y 


[3 • 6 • 3 • 6] 


2 


Y 


[3 • 122] 


4 


N 


[4 • 6 • 12] 


3 


Y 


[4 • 82] 


3 


Y 
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IV. A RIGOROUS LOWER BOUND ON W{k, q) FOR ARCHIMEDEAN A 



In Refs. we derived rigorous lower and upper bounds on Vr(A, q) for the triangular 

and honeycomb lattices, using a coloring matrix method of the type first applied by Biggs to 



obtain such bounds for the square lattice ||3J] . We showed that these upper and lower bounds 
rapidly approached each other for large q and became very restrictive even for moderately 
large q. Furthermore, we found that for a wide range of q values, the lower bounds were 
very close to the respective actual values of W{K,q). Accordingly, here we shall focus on 
rigorous lower bounds for W{K, q) for the general class of Archimedean lattices, of which the 
homopolygonal lattices are a special case. We shall derive a general lower bound applicable 
to any Archimedean lattice. In addition, we shall derive lower bounds for the duals of 
Archimedean lattices. 

Before proceeding, it is necessary to recall a subtlety in the definition of the function 
W{K,q). As we discussed in Ref. 0], the formal eq. ( [L.lD is not, in general, adequate to 



define Vr(A, q) because of a noncommutativity of limits 

lim lim P(A,g)^/" ^ lim lim P(A,g)^/" (4.1) 

n^oo q—>q3 q—*qs n^oo 

at certain special points g^- At these points, one must also specify the order of the limits in 
( [4.1|) . We denote these definitions as 



Wi{G},q.)n,,. ^ lim \im P{G,qY^- (4.2) 



and 



W{{G},q.)D„^ ^ lim lim P(G,g)i/" (4.3) 

where {G} denotes the n oo limit of the family of ri-vertex graphs of type G. One 
can maintain the analyticity of W{{G},q) at the special points of P{G,q) by choosing 
the order of limits in ([4. 2]); however, this definition produces a function iy({G}, g)_D,„ whose 
values at the points differ significantly from the values which one would get for P{G, qsY^^ 
with finite-n graphs G. The definition based on the opposite order of limits, ( [4.3| ) gives the 
expected results like W{{G},qs) = for q^ = 0,1, and, for G ^ A (i.e., for G containing 
at least one triangle as a subgraph), also qs = 2. However, this second definition yields a 
function iy({G}, q) with discontinuities at the set of points {qs}- Following Ref. Jl^, in our 
results below, in order to avoid having to write special formulas for the points qs, we shall 
adopt the definition Dqn but at appropriate places will take note of the noncommutativity 



of limits (|4.1|). As will be evident from the derivation, our rigorous lower bounds are on the 
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function W{A,q) = Vr(A,g)£)g„, and they apply for the range of (positive integer) q values 
for which the relevant coloring matrices (see below) are nontrivial, i.e., for q > xi-^)] the 
values of x(A) were listed in Tables H and |I|. 

Clearly, a general upper bound on a chromatic polynomial for an n-vertex graph G is 
P{Gi<l) < (f"- This yields the corresponding upper bound W{{G},q) < q. Hence, it is 
natural to define a reduced function that has a finite limit as g ^ oo, 

Wr{{G},q) = q-'W{{G},q) (4.4) 

When calculating large-g Taylor series expansions for W functions on regular lattices , it is 
most convenient to carry this out for the related function 

for which the large-g series can be written in the form 

oo 

W{K,y) = l+Y.w^,my"' (4.6) 

m=l 

with 

y = ^ (4.7) 

For duals of Archimedean lattices one can use the same formulas with the replacement 
A ^e//, where Ae// was given in eq. ( p.ll| ). 
Our rigorous lower bounds are of the form 

iy(A,g) >iy(A,g), (4.8) 

where the subscript £ denotes "lower", and we shall give the explicit expressions for W{K, q)i 
below. We shall use two other equivalent forms of the bounds, namely on the functions 
Wr{A, q) and TF(A, y), both of which have the finite limit Wr{k, q = oo) = W(A, y = 0) = 1 
and hence are convenient to compare with large-g (small-?/) series. Thus, for the latter 
function, the bounds read 

W{k,y)>W{k,y), (4.9) 
where the reduced lower bound function is defined by analogy with ( [4.5| ) as 

Our general rigorous lower bound, which is a major result of the present paper, is proved 
in the following theorem: 



14 



Theorem 



Let A = {YliPi') be an Archimedean lattice. Then for (integer) q > ^i-^^Q) = 

W{A, q)D„„ has the lower bound 



w{(npru),= ^-'';:f' (4.11) 



where the {i} in the product label the set of pj-gons involved in A, Up- was defined in eq. 
and 



s=0 



where P{Ck, q) = {q — 1)*^ + (— l)'^(g — 1) is the chromatic polynomial for a A;- vertex circuit 
graph, i.e., polygon. This lower bound takes a somewhat simpler form in terms of the related 
function W{K, y)f. 

w{{\{pT).y), = n[i + i-^ry'-'T' ^4.13) 



Proof 

We consider a sequence of finite 2D Archimedean lattices with periodic boundary condi- 
tions in one direction, say the x direction, and either periodic or free boundary conditions in 
the orthogonal, y, direction. Denote the lengths of the finite lattice in these two directions 
as m and n and the finite lattice of type A as A^xn- Extending the method that Biggs 
used for the square lattice [^, we introduce a coloring matrix T, somewhat analogous to 



the transfer matrix for statistical mechanical spin models. The construction of T begins by 
considering an n-vertex path £„ along a vertical path on A^xn- For the square lattice, it 
is obvious what is meant here, and one can also represent the other homopolygonal lattices 
as square lattices with bonds added or deleted (for the triangular and honeycomb lattices, 
respectively; see Fig. 1 of Ref. \TE\ and the discussion in Ref. [0]). Thus, as is well known. 



the triangular lattice can be deformed so that it is represented as a square lattice with bonds 
added to connect the lower left and upper right vertices of each square, and the honeycomb 
lattice can be deformed to make a brick lattice, with the long side of the bricks oriented 
vertically. These deformations do not affect the coloring properties of the lattices. With 
these representations of the triangular and honeycomb lattices, one can choose a vertical 
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path in the same manner as for the square lattice. Analogous paths can be defined on het- 
eropolygonal Archimedean lattices (an illustration is given in Fig. |l| for the (3^ ■ 4 ■ 3 ■ 4) 
lattice). In all cases except the (4 ■ 6 ■ 12) lattice, these are simple (i.e., unbranched) chains. 
For the (4 ■ 6 ■ 12) lattice, the paths are chains with single-bond branches. We shall give the 
proof first for the seven Archimedean lattices where the paths are simple chains which do not 
intersect, namely, the three homopolygonal lattices and the (3^-4^), (3^ • 4 • 3 ■ 4), (3 ■ 4 ■ 6 ■ 4), 
and (4 ■ 8^) heteropolygonal lattices. We then give the proof for the three remaining lattices 
where the paths do intersect, and, separately, for the special case of the (4 ■ 6 ■ 12), where 
the paths, although non-intersecting, are branched chains. The number of allowed colorings 
of the path is P{Cn, q) = q{q — 1)"^^ = Af. Now focus on two adjacent paths and C'^. 
Define compatible g-colorings of these paths as colorings such that no two vertices v & Cn 
and v' G C'^ connected by a bond of A^xn have the same color. One can then associate with 
this pair of paths an A/" x A/" dimensional symmetric matrix T with entries Tc„,c'„ = 1 or if 
the g-colorings of and C'^ are or are not compatible, respectively. Then for fixed lengths 
of the lattice in the x and y directions, m and n, P{A„ixn, q) = Tr{T^). For a given n, since 



T is a nonnegative matrix, one can apply the Perron-Frobenius theorem ||3^ to conclude 
that T has a real positive eigenvalue Xmax,n{q)- Hence, for fixed n, 

^im Tr(T'")i/(™")-.Al", (4.14) 

so that 

iy(A,g) = jim Al", (4.15) 

Denote the column sum Kj{T) = J^iLiTij (equal to the row sum Pj{T) = YjiLiTji since 
T'^ = T) and S{T) = Y.f,j=iTij] note that S{T)/U is the average row (column) sum. The 
lower bound is then a consequence of the (r = 1 special case of the) theorem that for a 
nonnegative symmetric matrix T [36| \max{T) > [S{T'^)/J\f]^^'^ for r = 1,2, ... together with 



eq. (|415|) : 

W{A,q) > W{A,q), = lim(^^y^ (4.16) 

This is the general method; we next proceed to calculate S{T). To do this we observe 
that the adjacent paths Cn and C'^ can be chosen such that the strip-like section of the 
lattice between them contains each of the types of polygons comprising A (this essentially 
amounts to the orientation of the lattice to define the vertical direction). Each polygon pi 
has the chromatic polynomial given above as P{Cp.,q). Next, we use a basic theorem from 



graph theory: if G and G' are graphs that intersect in a complete graph Kr then 
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P{G U G',q) = P{G,q)P{G',q)/P{Kr,q), where P{Kr,q) = IXZliq - s). We then apply 
this theorem, taking into account that the number of polygons of type pi per vertex is Up., 
to obtain the result that for paths of length n, S{T) = c{q) Yli Dp.{qY^i'^'^^^\ where h is an 
unimportant constant independent of n, and the prefactor c{q) = q{q — l){q — 2) if A contains 
triangles, A 3 A, and c{q) = q{q — 1) otherwise. Then taking the n —>■ oo limit in ( [4 .161) 



yields eq. (|4.11| ). Equivalently, using the definition of the reduced function ( [4 .51) , we obtain 



( [4.13| ). For the three lattices, viz., (3^-6), (3 ■ 6 ■ 3 ■ 6), and (3 ■ 12^), where the adjacent 
paths Cn and C'^ do intersect, the calculation involves a technical complication due to these 
intersections: the expression for S{T) involves an additional factor of {q — 1)^^, where A is 
the fraction of points on each path that coincide with points on but this factor is 
exactly cancelled by the same additional factor in the expression for A/" in the denominator of 
the ratio S'(T)/A/'. The additional factor in the denominator arises because one must correct 
for the undercounting of the points on the paths £„ due to the intersections. Because the 
factor cancels, the result is the same for as for the lattices with nonintersecting adjacent 
paths Cn and C'^. 

Finally, we consider the (4 ■ 6 ■ 12) lattice. Let us consider the lattice as oriented so that 
the 12-gons have vertical left and right bonds and horizontal top and bottom bonds, with a 
hexagon vertically above the 12-gon and adjacent to it, and label the vertices of each 12-gon 
in a clockwise manner starting with the left vertex on the top bond. The path that we 
use goes from vertex 6 of a given 12-gon to vertices 5, 4, 3, and then crosses over to vertex 
8 of the adjacent 12-gon in the upper right ("northeast") direction, has a 1-bond branch to 
vertex 7 of this 12-gon, and then continues along on vertices 8 through 12 and then 1, with 
a 1-bond branch to vertex 2, continuing from vertex 1 upward to vertex 6 of the vertically 
adjacent 12-gon in the northwest direction, and so forth. The strip enclosed by two adjacent 
paths of length n of this type has a chromatic polynomial 

S{T) = q{q - 1) [D,{q)'D,{qfD,2{q) + 2D^{q)D^{qYD^{q) + qD^^qfD^^q)]^'^'' (4.17) 

The chromatic polynomial for the path of length n is the same as that for an unbranched 
path, since both are tree graphs. Hence, for this case, the lower bound actually reads 
W{{A ■ 6 ■ 12), g) > W{{A ■ 6 ■ 12), g)^/, where 

,,,,,, , , [D,{qfDM^D,M + '^DMD,{qfD,{q) + qD,{qfD,{q)\^'' 

4 ■ 6 ■ 12 , g = ^ ^ 4.18 

q - 1 

For the above-mentioned range of q under consideration here, i.e., q > ■ 6 • 12)) = 2, the 
two additional terms in the square brackets are positive, so that 
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iy((4-6-12),g), 



q-1 



(4.19) 



which is of the form ( [4.11 ). (As we shall discuss below, the difference between the bounds 
( [4.18| ) and ([4.19|) is very small.) This completes the proof. 



For the homopolygonal Archimedean lattices A = (p") (with a = A given by (|2.5| ) and 



Up = 2/{p — 2) from (^.7D), our bound ( [4.11| ) reduces to 



W{ip^),q)^ 
or equivalently, ( |4.13| ) reduces to 



q-1 



.y 



1 + {~iru 



2/(p-2) 



(4.20) 



(4.21) 



In Table pT| , we list the explicit forms of the lower bounds Vr(A, y)£ for the Archimedean 
lattices, including a comparison with small-y series, to be discussed below. 

TABLE III. Rigorous lower bounds W{A,y)i for Archimedean lattices A = (YliPi'') given by eq. 

( 4.13| ). For homopolygonal Archimedean lattices A = (p^), the Taylor series expansion of IF(A, y)g 

coincides to order 0{y'"') with the series expansion of IF(A, y), where ic = imax = 2(p — 1), the 

values of which are listed in the third column, and the subscript c stands for "coinciding". For 

other lattices, the series expansions of W{A,y)£ and W{A,y) coincide to at least order 0{y^''). 



A 


W{K,y), 


(36) 




(44) 




(63) 


(l+y5)l/2 


(3^ • 6) 


(1- 2/2)4/3(1 +y5)l/6 


(33 . 42) 


(l-y2)(l + y3)l/2 


(32 • 4 • 3 • 4) 


(l-2/2)(l + y3)l/2 


(3 • 6 • 3 • 6) 


(l-y2)2/3(i+y5)l/3 


(3 • 4 • 6 • 4) 


(l_y2)l/3(i+y3)l/2(i^y5)l/6 
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(3-122) (l_y2)l/3(i+yll)l/6 13 

(4 -6 -12) (l + y3)i/4(i + y5)i/6(;L + yii)Vi2 11 

(4-82) (l + :V^)V^(l + ?y7)V^ 12 



19 



We give some illustrations of the method of the proof here. First, we illustrate this for a 
lattice where the paths do not intersect, namely, the (3^ ■ 4 ■ 3 ■ 4) lattice shown in Fig. |l[ 
The paths and CJ^ are depicted by the thicker lines, and the strip between these consists 
of a sequence of squares and double triangles. For technical convenience, we consider a path 
to start at the lower lefthand vertex of a square, and we take n to be odd, so that the path 
length is even. The sum S'(T) is calculated starting from the basic graph G433 comprised of a 
square and two adjacent triangles (say above the square) for which the chromatic polynomial 
is -P(G433, g) = g(g — 1)1^3(5)^1^4(0'). For a strip lying between the adjacent paths and 
CJ^ (both of length n — 1), starting the count from the lower left corner of a square, there 
are r = (n — l)/2 6*433 graphs, so, in an obvious notation, 

S[T) = P{Gl,,,q) = q{q - 1) [/^3(g)^ ^(g)]^"-^)/^ (4.22) 

Dividing by A/" = q{q — 1)"^^ and taking the n — 00 limit of the 1/n'th power of the ratio 
as in eq. (|4.16| ), one obtains the rigorous lower bound 



W^((3^.4.3.4),g), = ^^^M^^iM^ (4.23) 

q- 1 



which is seen to be the special case of (|4.11| ) for this lattice. Another explicit example is 



provided by the (4 ■ 8^) lattice, which we have previously discussed . 



To illustrate the proof in a case where the paths £„ and £^ have a nonzero overlap, we 
discuss the (3 • 6 ■ 3 • 6) lattice, shown in Fig. |^. To visualize the path start at the vertex 
A at which two triangles touch and move "northwest" to the leftmost vertex of the hexagon 
above these triangles; then move northeast to the vertex common to the two triangles above 
the hexagon, then northwest again to the leftmost vertex of the next higher hexagon, etc., 
continuing upward in this zigzag manner. For C'^, start from the same vertex, but move 
first northeast to the rightmost vertex of the hexagon above the triangles, then northwest 
to the vertex common to the two triangles above the hexagon, and so forth. For technical 
convenience, we consider the paths to be of length a multiple of 4. Evidently, the n-vertex 
paths Cn and C'^ so defined have a nonzero overlap (intersection) consisting of the vertices 
that are shared in common by the pairs of triangles traversed along the route. For consistency, 
we assign alternate intersection points to be on adjacent paths. These intersection points 
comprise 1/3 of the total number n of the vertices along each path. The strip between the 
paths Cn and C'^ consists of a sequence of triangle-hexagon-triangle graphs, which we denote 
as 6*363. For a path starting from a vertex of type A, there are r = {n — l)/3 linked G363 
graphs in this strip. The chromatic polynomial of the strip is 

SiT) = P{Gl,,, q) = q[{q - l)D,{qf D.iq)^-'^'' (4.24) 
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However, in evaluating the denominator for eq. ([4.16|) , it is necessary to correct for the 
vertices that are common to both paths; taking into account that these comprise a third 
of the total points on each path and including this correction factor yields the denominator 
M = q{q — Evaluating eq. (|4.16| ) then gives the resulting lower bound 



V^((3.6.3.6),^^5(«^ ,4,25) 



which again is the special case of ( |4.11| ) for this lattice. 



Finally, we discuss the difference between the lower bounds (|4.18| ) and ( [4.19|) for the 



(4 ■ 6 ■ 12) lattice. This difference shows that, at least for this lattice, the special case ( [4.19 ) 



is not an optimal lower bound. However, the numerical difference is extremely small. At 
q = x((4 ■ 6 ■ 12)) = 2, the bounds (|4.18|) and ([4.19|) are identical and equal to 1, which 
is also the exact value of W^((4 ■ 6 ■ 12), g = 2)D„q (see eq. ([4.33|) below). At g = 3, the 
difference between the bounds ( [4.18[ ) and ( [4. 19] ) is 1.3 x 10"^, and this difference decreases 
monotonically for larger q. This reflects the fact that the large-g series expansions of the 
corresponding reduced lower bound functions are identical to 0{q~^^): 

Wr{{A •6-12), q)e - Wr{{4 ■ 6 ■ 12), g), = ^g-^^ + 0{q-'') (4.26) 

We comment on some general properties of our rigorous lower bound on Vr(A, g) in 
eq. ( [4.11|) , and the consequent bound on the related functions Wr{A,q) and W{A,y) for 



Archimedean lattices A. First, 

w{il[pT),q), ~ ~ g, as g^oo (4.27) 

i g 

where X^i ^^iiPi — 2) = 2 as a consequence of eq. (p^). Hence, limg^oo Wr(^{IliPi ')} o) ^ = ^■ 

Second, note that when we compare our lower bound on Wr{A,q) to the large-g Tay- 
lor series expansion on each lattice A, this comparison is obviously restricted to the range 



of g for which these series expansions are applicable. As we discussed in Ref. [|^, for 
a given graph G, when one generalizes g to a complex variable, the reduced function 
Wr{{G},q) = q~^W{{G},q) is analytic in the g plane except on the union of boundaries 
B. These boundaries can separate the g plane into various regions Ri such that one cannot 
analytically continue Vr({G}, g) from one region to another. We defined the region contain- 
ing the positive real g axis from a minimal value gc(A) to g = oo as Clearly, the large-g 
Taylor series for iyj.({G}, g) apply in this region. Here, qdG) is thus defined as the maximal 
(finite) real value of g where W{{G},q) is nonanalytic. All of these considerations apply 
for the case of regular lattice graphs A defined, as above, as the n — > cx) limit of finite A„ 
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lattices. As our general study and explicit exact solutions in Ref. showed for g G -Ri, the 
limits ( [4.1|) commute and hence the two definitions (4^) and ([4.3|) coincide. Hence, in our 
discussions of series, we shall not need to distinguish between these two different definitions 
of W{A,q). Furthermore, in this region Ri, just as one can extend the definition of the 
function 1V(A, q) from integer to real q, so also one can carry out the same extension of the 
lower bound W{A, q)e. Hence, for q > gc(A), the lower bound (|4.11| ) applies for real, and not 
just integer q. In general. 



X(A) < g,(A) 



(4.28) 



In Refs. 161 and we showed that qdi^^)) = 4, gc((4^)) = 3, and gc((6^)) = (3 + V5)/2. 

Third, W{A, q)i is a monotonically increasing function of real q for q > x(A) |]38| . 

Fourth, if two different lattices involve the same set of polygons and have the property 
that each type of polygon pi occurs an equal total number of times as one makes the circuit 
around each vertex, then our lower bound ( [4.11D is the same for both. This is true of the 
(3^ ■ 4^) and (3^ ■ 4 ■ 3 ■ 4) lattices, for which the lower bound is given by ( [4.23|) . 

Fifth, it should be remarked that one can often choose different types of paths and 
C^. For example, in Ref. |18| we derived a lower bound W{{6^), q)e^dimer = iQ~ l)'^^"^ / q^^"^ by 
choosing paths consisting of certain dimers. In Ref. [O] , using paths of the type described 



here, we obtained the lower bound of the form (|4.11| ) 



W{i6'),q)t 



q-1 



> W q)i,dimer 



(4.29) 



which, as indicated, is slightly more stringent. 

Finally, the fact that for the circuit with an even number, say 2k, of vertices, the chromatic 
polynomial is P{C2k, ? = 2) = 2, together with the definition ( 4.12 ), it follows that |39 



D^kiq = 2) = 1 



(4.30) 



Hence, for bipartite Archimedean lattices A^jp , where our lower bound ( [4.11D can be applied 

at g = xi-^bip.) = 2, we have 



W{Abip.,q = 2) 



(4.31) 



Thus, in this case, although our bound refers to W{A,q)Dq„, it has the same value as 
W{Af,ip,, q = 2)d„,, as given below in eq. ( [4.33| ). 

Note that for lattices Auqc that are uniquely g-chromatic (UQC) for q = x(^)) it follows, 
a fortiori, that 
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W{KuQC.q = xmD^, = l (4.32) 
Thus, in particular, for a bipartite lattice 

W{k,,p,,q = 2)o^^ = l (4.33) 
and for the tripartite lattices (3^) = [6^], (3^ -6), [4 ■ 8^], and [4 ■ 6 ■ 12], 

iy(At„p.,g = 3)z5„, = 1. (4.34) 

V. A RIGOROUS LOWER BOUND ON W{K, q) FOR DUAL ARCHIMEDEAN A 

Using methods similar to those that we used for the Archimedean lattices, we have also 
calculated rigorous lower bounds for the duals of these lattices. As with the lower bounds for 
Archimedean lattices, our bounds for the dual lattices again apply for q > x(A) for each A. 
For seven of the lattices, we are able to obtain the respective lower bounds by using paths 
Cn that either do not intersect, in the case of the homovertitial lattices and also for [3^ ■ 6], 
[3^ -42], [32 ■ 4 ■ 3 ■ 4], and [3 ■ 4 ■ 6 ■ 4], or intersect, in the case of the [3 • 6 • 3 • 6] lattice. For 
these, by the same methods as before, we obtain the general rigorous lower bound 

wm^h<i)i = ^^ (5.1) 



where p and Up were given in eqs. (|2.9| ) and ( p.lO| ). Using eq. ([4.10|) with A replaced by 
Ae// given in eq. ( ^.llj ), this is equivalent to 

W{[Y[v':%Q)i={^ + i-W-T (5.2) 

i 

The proof is the same as the one given before for Archimedean lattices, with the simplification 
that only one kind of polygon is involved in the calculation of S{T). The homovertitial 
Laves lattices have already been dealt with as homopolygonal Archimedean lattices; for the 
heterovertitial Laves lattices satisfying the above condition on paths we find the specific 
bounds 

Wi[3' ■ 6], q)e = Wi[3' ■ A% q), = Wi[3' ■ A ■ 3 ■ A], q)e = (5.3) 

q - 1 



and 



W{[3-6-3-6],q), = W{[3-4.6-4],q),= ^=^^^-^^ (5.4) 
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M.,|4.8^1.,)«^i^^(i^ (5.5) 

q — 1 Q — i- 



We note that the lower bound ( |5.4| ) is the same as the one that we derived for the square 
lattice (4'^) = [4^], and the lower bound ( p.5| ) is the same as the one that we derived for the 
triangular lattice, (3^) = [6^]. 

However, for one of the lattices, viz., [4 ■ 8^], where such paths exist, we have been able 
to obtain a more stringent lower bound by using a different kind of coloring matrix method, 
in which this matrix is defined in terms of the compatibility not of adjacent paths, but 
of adjacent chains of graphs that are more complicated than paths. Consider the [4 ■ 8^] 
lattice shown in Fig. |^. For this lattice, one obtains the lower bound ( p.5| ). But rather than 
paths, one can consider chains of graphs constructed as follows: consider a vertical line on 
the lattice, and let denote the chain of triangles whose bases include vertices from 1 
to n along this vertical line and whose apexes point to the right. Evidently, each pair of 
sequential triangles intersect at the common vertex on their bases. Define Ai'^ to be the 
corresponding equal-length chain of triangles whose bases form the unit segments along the 
next (adjacent) vertical line of the lattice to the immediate right of the previous one. For 
technical convenience, let n be odd. Then define a new type of (symmetric) coloring matrix, 
of dimension Af x Af, where Af = P{M.n,<i) = P{.-M'n,q), with entries Tm^m'u = 1 or if 
the g-colorings of the chains of graphs M.n and AA'^ are or are not compatible, respectively. 
One then proceeds as before, to derive the lower bound ([4.16|) . In the present case, 

71 — 1 

P(^„,g)=g[(g-l)(g-2)r , r = ^— (5.6) 



and 



77 — 1 

5(T) = g(g-l)[(g-2)2(g2-5g + 7)r , r=^ (5.7) 



so that, evaluating the limit in ( [4.16| ), we obtain 



(g-2)(g2_5g + 7) 



1/2 



(5.^ 



For the relevant range, g > ■ 8^]) = 3, the lower bound (|5.8|) always lies above (|5.5| ), i.e., 
is more stringent. We list the corresponding lower bound iy([4 ■ 8^], y)^ in Table IV). 

Using this different coloring matrix method based on chains of graphs that are not sim- 
ple paths, we similarly obtain the following rigorous lower bounds for the other two dual 
Archimedean lattices: 

r n 2/3 

(g-2)(g-3) 

W{[?> ■ 12% q), = ^q \y/3 (5-9) 
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„.(|4.,.1.].,),.(-^|1Z^ (5.0, 

The corresponding lower bounds W^(A, y)i are listed in Table To retain the symmetry 
between the Archimedean and Laves lattices, we also list the lower bounds for the square, 
triangular, and honeycomb lattice in their respective Laves forms, [4"^], [6^], and [3^]. 

TABLE IV. Rigorous lower bounds W{A,y)g for duals of Archimedean lattices, A = [Hj^i']) 
where Vi = Aj. For homovertitial duals, A = [v^] = [A^] where p and A are related by eq. 
(|2.6| ), the Taylor series expansion of W{A,y)£ coincides to order 0{y^'') with the series expansion 
of 1^(A, y), where ic = imax = 2(p — 1), the values of which are listed in the third column, and 
the subscript c stands for "coinciding". For other lattices, the series expansions of W{A,y)i and 
W{A,y) coincide to at least order Oiy^"). 



A 


W{A,y), 


[A 


(l+y5)l/2 




l + y3 


[63] 




[3^ • 6] 


(1 - y'?" 


[33 . 42] 


(1 - y'?" 


[32 • 4 • 3 • 4] 


(1 - y'?" 


[3 • 6 • 3 • 6] 




[3 • 4 • 6 • 4] 




[3 • 122] 


(l + 2/)2(l- 32/ + 22/2)2/3 


[4 • 6 • 12] 


(l+2/)2(l-2/)(l- 32/ + 32/2)1/3 


[4 • 82] 


(l + 2/)2(l-2/)i/2(i_3y + 3y2)i/2 



25 



VI. LARGE-g SERIES 



A priori, a lower bound on a given function need not agree to any, let alone many, 
terms in the Taylor series expansion of that function about some special point. In our earlier 
comparison of rigorous lower bounds on W{A, y) with large-g (i.e., small-?/) series expansions 
for the respective exact W{K,y) functions for the homopolygonal lattices A = (p^), where 
A = 2p/(p — 2), we found [p!7| -p!9| that our lower bounds and the series coincided to order 
O (?/*'=) where 

ic = imax = 2(p - 1) (6.1) 

In the case of the (6^) (honeycomb) lattice, our lower bound thus coincided with the small-y 
series for the exact function to 0{y^^), i.e., the first eleven terms. Thus, in addition to being 
a rigorous lower bound, our expression for W{{6'^),y) became an extremely accurate approx- 
imation to the exact function, W{{6'^),y) for even moderately large q. We also considered 
one heteropolygonal Archimedean lattice, the (4 ■ 8^) lattice, and calculated both a lower 
bound and small-y series for PF((4 ■ 8^), y) ||T^; in this case we found that our lower bound 
coincided with the first thirteen terms of our series. 

Here we report our calculations of small-?/ series for W^(A, y) on the remaining seven 
Archimedean lattices and compare these with our rigorous lower bounds, thereby achieving 
a complete comparison for all Archimedean lattices. Our calculations use the methods of 
Ref. Since our main purpose was an exploratory study of the extent to which our lower 
bounds coincided with the series expansions, we have not attempted to compute the latter 



to very high order; it is straightforward to carry these expansions to higher order (e.g. ||41|| ). 
Our results are listed below: 

W{{3' ■ 6), ?/) = 1 - ^y' + \y' + 0(?/^) (6.2) 
TF((33 . e),y) = W{{3' . 4 ■ 3 ■ 4), ?/) = 1 - ?/2 + ^ + . + 0{y') (6.3) 

2 3 4 

TF((3-4-6-4),?/) = l-|- + |--|- + 0-?/^ + O(?/^) (6.4) 
W((3 ■ 6 ■ 3 ■ 6), ?/) = 1 - ^?/2 - + ^y' - - ^y' - + 0(?/9) (6.5) 

W((3 ■ i2%y) = i-\y'- \y' - ^,y' - fy - f^' + - - + 0{y'') 

(6.6) 
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W((4 . 6 ■ 12), y) = 1 + + 1 ^5 _ + ^^8 + _7_^9 _ ^^10 + + ^(y^) (6.7) 

W((4 ■ 8^), y) = 1 + - ^y' + \i + ^y' + - + 0{y'') (6.8) 

We have also calculated low-order series for dual Archimedean lattices in order to make 
an initial comparison with our rigorous lower bounds. We have compared each of these 
series with the corresponding small-?/ Taylor series expansions of our rigorous lower bounds 
in Tables pT| and fV[ We find that the latter expansions coincide at least to the order Oly^"), 
where the respective values of ic are listed in Table |T| for each lattice. As is evident from 



this table, the striking fact that we showed previously |T8| , p!9| , viz., that the series expansions 
for W{A,y)£ coincide to very high order with the respective series expansions of the exact 
functions W{A, y) for A = (6^) (to 0(y^°)) and for A = (4 ■ 8^) (to at least O^y^"^)), is not 
restricted to just these lattices. Indeed, we observe that for both of the lattices (3 ■ 12^) and 
(4 ■ 6 ■ 12), the respective small-y series for the lower bound and the exact function coincide 
to at least 0{y^^) and 0{y^^), respectively Thus, in general, this establishes that for a 
number of lattices our rigorous lower bounds actually serve as quite good approximations to 
the exact W functions, especially for large q. 



VII. PLOTS OF WriKq) AND W{k,y) 

In Fig. 1^ we plot Wr{A,q) = q~^W{A,q) for the eleven Archimedean lattices, for y = 
l/(g — 1) in the range < y < 0.30, corresponding to q greater than about 4. We have 
evaluated these from our small-y series expansions of W^(A, y) together with the definition 
([4.5|). We find that Wr{A, q) and W{A, q), and hence the ground state entropy of the g-state 
Potts antiferromagnet, S'o(A,g) = kBhiW{A,q), are monotonically decreasing functions of 
the coordination number A(A) of the lattice A. This is a consequence of the fact that as one 
increases A, one increases the number of constraints restricting the coloring of each vertex of 
the lattice. Indeed, one can observe, especially for small y, that the eleven curves in Fig. || 
fall into four groups, for the four values of A = 3, 4, 5 and 6. For further analytical purposes, 
we include also a plot of the reduced function W{A, y) in Fig. ^ Again, this is calculated 
from our small-?/ series. From our detailed comparison of lower bounds and small-y series for 
the homopolygonal lattices and the (4 ■ 8^) lattice with Monte Carlo calculations, we expect 
that these curves are accurate over the range shown, i.e., y = to y = 0.3. We have also 
checked this by evaluating the sizes of the highest calculated terms in our series. 
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FIG. 6. Plots of Wr{A,y) = q ^W{A,y) as a function of y = — 1) for Archimedean lattices 
A. For visual clarity, the curves are shown alternately as solid (s) and dashed (d). The order of 
the curves, from bottom to top, is: (3^) (s), (3^^ • 6) (d), (3^ • 4p) (s), (3 • 6 • 3 • 6) (d), (3 • 4 • 6 • 4) (s), 
(4^) (d), (3 • 12^) (s), (6^) (d), and (4 • 8^) and (4 • 6 • 12) as a single solid curve to this resolution. 
See text for further details. 
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FIG. 7. Plots of W{A,y) as a function of y = l/(g — 1) for Archimedean lattices A. For visual 
clarity, the curves are shown alternately as solid (s) and dotted (d). The order of the curves, from 
bottom to top, is: (3^) (s), (3^ • 6) (d), (3^ • 4^) (s), (3 • 6 • 3 • 6) (d), (3 • 122) ^^^^ (3 . 4 • 6 • 4) (d), 
(6^) (s), (4 • 8^) and (4 • 6 • 12) as almost coincident dotted and solid curves, and finally, (4*^) (d). 
See text for further details. 

We note the following theorem: 

Theorem: If a lattice A' can be obtained from another, A, by connecting disjoint vertices of 
A with bonds, then W{A', q) < W{A, q) (where g > is an integer). For 
q < max{gc(A), gc(A')}, this inequality applies to the W functions defined with the Dng 
order of limits in eq. ( [4 .31) , while for larger q, it applies to either of the definitions ( [4 .31) and 
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( [4.2|) since they are equivalent for this latter range of q. 



Proof. Consider a finite, n-vertex lattice of type A, denoted A„. We use the 

addition-contraction theorem from graph theory, the statement of which is the following: 

let G be a graph, and consider any two vertices f , v' on G which are not adjacent, i.e., not 

connected by a bond of G. Denote the graph with a bond connecting v and v' as 

and the graph with these two vertices identified as Gi,=„'. Then 

P{G, q) = P{G,j^b^^',q) + P{G^=^', q). Since P(G^=^/, q) > 0, it follows that 

P{Gv-b-v', q) < P{G^ q)- Now take A„ = G and add bonds as necessary to construct the 

lattice A^. Each time one adds a bond, one gets an inequality on the correponding 

chromatic polynomials, thereby producing a sequence of such inequalities, expressing the 

fact that the chromatic polynomial for the original lattice is > that for the lattice with one 

bond added, which is > that for the lattice with two bonds added, etc. Together, these 

yield the inequahty P(A„, q) > P{A'^, q). Now let n ^ oo to obtain 

W{K^q)j:,^^ > W{A',q)£)^^. For q > max{gc(A), gc(A')}, both of these definitions are 

equivalent, so one can drop the subscripts D^q and Dg^. 

We give two examples of this theorem. (Subscripts indicating orders of limits in the 
definition of W are understood where necessary.) First, as recalled above, the square lat- 
tice can be obtained from the honeycomb lattice by such bond addition, and consequently, 
W{{4:^),q) < W{{6^),q). Second, the triangular lattice can be obtained from the square 
lattice by bond addition, so W{{3^),q) < W{{A^),q). 

Concerning the plot of W{A,y), we observe one basic feature: as y increases from 0, i.e. 
q decreases from oo, the initial behavior of W{A, y) can be understood as resulting from the 
leading term in the small-?/ expansion 



where g denotes the girth of A. That is, for even (odd) g, W{K,y) increases (decreases). 

It is of interest to compare our results for the dependence of the W^(A, q) coloring function, 
or equivalently, the exponent of the ground state entropy of the Potts antiferromagnet, on the 
lattice coordination number A with the A-dependence of other models that exhibit ground 
state entropy. We first consider the (spin 1/2) Ising antiferromagnet (lAF) on the triangular 
and kagome lattices. The exact solutions in these cases yield 



Ty(A,y) = l + (-l)%y^-i + ... 



(7.1) 



1 /•2- d?9 



ln(3 + 2P) = - / duj ln(2 cos uj) ~ 0.3231 (7.2) 
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for the triangular lattice, where 

P = cos(^i) + C0S(^2) + 008(^1 + 62) (7.3) 



and, for the kagome lattice |25 



SoilAF, (3 ■ 6 ■ 3 • Q))/kB = ^ j^^ ln(21 - 4P) ~ 0.5018 (7.4) 

Thus, although the ground state entropy is accompanied by frustration in these cases, in 
contrast to the g-state Potts antiferromagnet for the range q > x(^) considered in the rest 
of this paper, it is again true that the ground state entropy decreases as the coordination 
number A of the lattice increases. A similar dependence has been reported in the case of the 
quantum Heisenberg antiferromagnet; for the cases where this model involves frustration, 
recent studies indicate that it has a ground state with long range order (albeit nonmaximal) 
on the triangular lattice, but with nonzero ground state entropy and no long range order in 
the case of the kagome lattice . 

We also make a comparison with the A dependence of generalized ice models. At normal 
pressures p ~ 1 atm, ice forms a wurzite crystal with fixed coordination number A = 4, so, 
of course, one cannot vary A in a realistic ice model. A very accurate estimate of the entropy 
of ice was obtained by Pauling 0: So{ice)p^moiar/R = ln(3/2) = 0.4055, where, as above, 
R = NAvog.kB, to be compared with the measured value of Sq/R = 0.41 ±0.03 [|l|. Although 
one cannot vary A for real ice itself at p = 1 atm., one can consider an abstract statistical 
model constructed to have two-valued variables (say arrows) assigned to links subject to the 
constraint of local arrow conservation, inspired by the physical and chemical constraint of 
local electric neutrality in real ice. Clearly, such models can only be defined on a lattice 
with even coordination number A. One sees that this immediately constitutes a difference 
with spin models, which can be defined on lattices with even or odd A, as well as lattices 
such as the dual Archimedean lattices on which different vertices Vi have different degrees 
Aj. We recall that there is a straightforward generalization of the Pauling estimate for the 
exponential of the entropy of real ice to that for an abstractly defined ice-type model; this 
consists of the product 2^ describing the unconstrained number of positions of the hydrogen 
ions on the bonds (where E = {A/2)V and V denotes the number of vertices on the lattice), 
multiplied by a reduction factor which is the fraction of the number of configurations for each 
vertex (oxygen location) that satisfy local electric neutrality. For each vertex this fraction 
is (^^2)/^^' combining these two factors yields the generalized Pauling estimate for the 
exponent of the entropy, per site, for ice-type models (denoted by a subscript /), 

W,{A)p = 2-^/^(^\ (7.5) 
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This is a monotonically increasing function of (even) A. Since ( [7 .51) is known to be a rigorous 
lower bound on Wi{A) |^|, it follows that Wi{A) is also a monotonically increasing function 
of (even) A. This shows that different models that both exhibit nonzero ground state 
entropy, such as (i) the zero-temperature Potts antiferromagnets for q > x(^) various 
lattices A; the (frustrated) Ising antiferromagnet; and the (frustrated) quantum Heisenberg 
antiferromagnet on the triangular and kagome lattices on the one hand, and (ii) ice-type 
models on the other hand, can have quite different dependences on lattice properties such 
as the coordination number A. 



VIII. CHROMATIC ZEROS AND SUPPORT FOR A CONJECTURE FOR 

REGULAR LATTICES 

We have calculated chromatic polynomials for finite sections of heteropolygonal 
Archimedean lattices and for several Laves lattices, and in each case have calculated the 
zeros of these polynomials, i.e., the respective chromatic zeros for these lattices. We have 
also done this for the simple cubic lattice. Our main motivation for these calculations is to 
check whether the results are consistent with the conjecture that we made previously |P0|| , 



namely that a sufficient condition that Wr{{G},q) be analytic at 1/g = 0, (i.e. that the re- 
gion boundary B separating different regions where iy({G}, q) is analytic does not have any 
components that extend to complex infinity in the q plane) is that {G} is a regular lattice. 
(This is not a necessary condition; as our exact results in Ref. showed, there are many 



families of graphs that are not regular lattices but do have compact region boundaries B.) 
In Refs. ||T7| , pO| we compared the chromatic zeros for various families of graphs with exact 
results on the n — > oo limits of these graphs and the corresponding W{{G},q) functions, 
where {G} denotes the n oo limit of n-vertex graphs of type G. As we discussed, as 
n — s> oo, the chromatic zeros (aside from a well-understood discrete subset including zeros at 
g = 0, 1 and, for graphs containing at least one triangle, also q = 2) merge to form boundary 
curves B that separate regions of the complex q plane in which VTdG}, q) takes on different 
analytic forms. For finite graphs, these chromatic zeros lie near to, or, for some families of 



graphs [|T7| exactly on, the asymptotic boundary curves B. However, for cases where B does 
have components extending to complex infinity in the q plane, we found that the chromatic 
zeros calculated on finite lattices deviate strongly from the parts of B that extend to infinity 
(see also [0). In particular, as one increases the lattice size, one finds that the (complex- 
conjugate) complex chromatic zeros farthest from the real axis in the q plane move farther 
away from this axis. This, then, serves as a means by which one can infer, from calculations 
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of chromatic zeros on finite lattices, whether the corresponding region boundaries B have 
components that extend to complex infinity in cases where one does not have exact solutions 
for g) in n — s> oo limit available. As we have discussed before |ll7| , p0| , exact results 



for the triangular lattice and chromatic zeros calculated on finite square and honeycomb 
lattices fT^ yield boundary curves B that satisfy our conjecture. Summarizing our results 
for various Archimedean and dual Archimedean lattices, and for the simple cubic lattice, we 



find in all cases that the chromatic zeros are consistent with our conjecture in Ref. |2^. We 
show two typical examples in Figs. || and ^j. In both cases, we use free boundary conditions 
and choose sections of the lattice that have comparable lengths in the x and y directions. 
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FIG. 8. Chromatic zeros for a section of the (3 • 12^) lattice, with n = 48 vertices. 
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FIG. 9. Chromatic zeros for a section of the (4 • 8^) lattice, with n = 36 vertices. 

From our earlier comparisons of chromatic zeros for various families of graphs with the 
exact region boundaries B calculated in the limit of infinitely many vertices we know 

that it is possible to make some reliable inferences about these boundaries from the positions 
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of the chromatic zeros for finite lattices. Indeed, a subset of these chromatic zeros merge to 
form the boundaries B in this hmit. (There are also discrete isolated chromatic zeros such as 
those at g = 0, 1 and, if A contains triangles, also at g = 2.) Our calculations of chromatic 
zeros for Archimedean lattices and their duals in 2D and for the simple cubic lattice are 
consistent with the inference that in the thermodynamic limit, the respective boundaries 
(1) separate the complex q plane into at least two regions, one of which (denoted Ri in 
Refs. |1T7| , |20[| ) includes the positive real q axis extending to g = cxd and the circle at complex 
infinity, i.e., the image under inversion of the origin in the 1/q plane; and (2) the outermost 
component of the boundary B intersects the real q axis on the left at g = (for all A) and on 
the right at the point that we have denoted gc(A). Further studies on larger lattices will help 
to elucidate the detailed shapes of the boundaries B for various lattices. For example, using 
sufficiently large lattices together with comparisons of chromatic zeros for different lattices 
sizes to measure finite-size shifts of these zeros, one can carry out an extrapolation to the 
thermodynamic limit to determine the value of gc(A) with reasonable accuracy. Work on 
this is in progress. 



IX. CONCLUSIONS 

The subject of nonzero ground state entropy is a fundamental one in statistical mechanics. 
In this paper we have proved a general rigorous lower bound for W{A,q), the exponent of 
the ground state entropy of the g-state Potts antiferromagnet on an arbitrary Archimedean 
lattice. The function Vr(A, g) is also of considerable interest in mathematics, in particular, 
the coloring of (finite) graphs and their infinite-n limits. From calculations of large-g series 
expansions for the exact Vr(A, y) functions and comparison with our lower bounds on the 
various Archimedean lattices A, we have shown that the lower bounds are actually very good 
approximations to the exact functions for large g. We have also calculated lower bounds and 
series for the duals of Archimedean lattices. Finally, from calculations of chromatic zeros on 
a number of lattices, we have obtained further evidence for the conjecture that a sufficient 
condition for q~^W{A, g) to be analytic at 1/g = is that A is a regular lattice. 

This research was supported in part by the NSF grant PHY-93-09888. 
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